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Abstract 



I Recent experimental results on the color-suppressed nonleptonic decays — 

D^*^^TT^ provide evidence for a failure of the naive factorization model and for 
Oh. sizeable relative strong-interaction phases between class-1 and class-2 B — > D^*^7r 

(— I I decay amplitudes. The allowed regions for the corresponding ratios of (complex) 

isospin amplitudes and 02/ai parameters are determined. The results are inter- 
. preted in the context of QCD factorization for the related class-1 amplitudes in 

rS ' the heavy-quark limit. 



The problem of understanding nonperturbative strong-interaction effects in exclusive 
nonleptonic weak decays of liadrons has always been a challenge to theorists. Only in 
a few cases model-independent results based on controlled expansions in QCD can be 
obtained. In the absence of a quantitative theoretical description various attempts have 
been made to obtain simple, predictive parameterizations of decay amplitudes based on 
simple phenomenological assumptions. 

The most common of these approaches is the "naive" (or "generalized" ) factorization 
model, in which the decay amplitudes are estimated by replacing hadronic matrix ele- 
ments of four-quark operators in the effective weak Hamiltonian by products of current 
matrix elements determined in terms of meson decay constants and semileptonic form 
factors. "Nonfactorizable" strong-interaction effects are parameterized by phenomeno- 
logical coefficients a,, which depend on the color and Dirac structure of the operators but 
otherwise are postulated to be universal |, ^. One distinguishes class-1 and class-2 
decay topologies, which refer to the cases where a charged (class-1) or a neutral (class-2) 
final-state meson can be produced from the quarks contained in the four-quark opera- 
tors of the effective Hamiltonian. For instance, in decays based on the quark transition 
b —>■ cud mesons with quark content [ud) or {uc) can be produced in that way. The 
decay —>■ D~^tt~ is a class-1 process, in which the pion can be generated at the weak 
vertex, whereas — > D^ir^ is a class-2 process, in which the meson can be directly 
produced. The corresponding amplitudes are then expressed as 

A{B' ^ D+vr-) = z ^ V^,V:, {ml - ml) U F,^-^{ml) a,{D^) , 

V2A{B^ - D\^) = I ^ V,,V:, {ml - ml) fn F,^--{ml) a,{Dn) , (1) 

where F^^'^^q^) are B M form factors at momentum transfer g^. In other processes 
such as B^ D^TT^ both topologies can contribute and interfere. (Such processes are 
sometimes called class-3 decays.) In fact, isospin symmetry implies that 

A{B- D\-) = A{B^ D+TI-) + V2^(5° ^ DV) . (2) 

The large- iVc counting rules of QCD show that ai{DTx) = 0(1) and a2{Dn) = 0{1/Nc), 
which is why the class-2 decays are often referred to as "color suppressed" . 

In the naive factorization model one postulates that for a large class of energetic, two- 
body (or quasi two-body) B decays the coefficients ai and a2 are process-independent 
phenomenological parameters. These parameters are assumed to be real, ignoring the 
possibility of relative strong- interaction phases between class-1 and class-2 amplitudes. 
Surprisingly, despite their crudeness these assumptions seemed to be supported by ex- 
perimental data 0, ^. Within errors, the class-1 decays 5° D^*^+M- with M = 
n, p,ai, Ds, D* can be described using a universal value |ai| ^ 1.1 ± 0.1, whereas the 
class-2 decays B — > K^*^M with M = J/%jj,ijj{2S) suggest a nearly universal value 
\a2\ ~ 0.2-0.3 (which is more uncertain due to the uncertainty in the B — ^ K^*^ form 
factors). Moreover, the class-3 decays B^ D^*^^M~ with M = Tr,p, which are sensi- 
tive to the interference of the two decay topologies, could be explained by a real, positive 
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ratio 02/01 ~ 0.2-0.3, which seemed to agree with the determinations of |ai| and \a2\ 
from other decays. The missing hnk in this hne of argument was a direct measurement 
of 1 02 1 in the related class-2 decays 5° ^ L)(*)°M°. 

Recently, the idea of factorization in the class-1 decays 5° D^*^^L^, where L is 
a light meson, was put on a more rigorous footing. It was shown that the correspond- 
ing decay amplitudes can be systematically calculated in QCD in the limit where the 
decaying b quark is considered a heavy quark @, ||, ^ (see 0, H, P] for related earlier 
work). To leading order in A/mf, (with A a typical hadronic scale), but to all orders of 
perturbation theory, nonfactorizable strong-interaction effects can be described in terms 
of convolutions of hard-scattering kernels with the leading-twist light-cone distribution 
amplitude of the light meson L. The resulting QCD factorization formula allows us 
to compute the magnitude and phase of the parameters ai{D^*^L) systematically up to 
power corrections in Kjm^,. The values of a\ in different class-1 decays are not univer- 
sal, but the process-dependent corrections turn out to be numerically small [Q. Several 
types of power corrections to the ax parameters have been estimated and found to be 
small |Ty, [TI|, |r^. Hence, for the cases where L is a light meson there is now a solid 
theoretical understanding of the near-universal value |ai| ^ 1.1 observed experimentally. 

On the other hand, if the charm quark is treated as a heavy quark, then the QCD 
factorization formula does not apply for the class-2 decays & — > D^*^^L'^, and so the 
magnitude and phase of the a2{D^*^L) parameters are not calculable. The only nontrivial 
prediction in this case is that the class-2 amplitudes are power suppressed with respect 
to the corresponding class-1 amplitudes 0. The apparent universality of the \a2\ values 
extracted from experiment, and the absence of sizeable relative strong-interaction phases 
between the various Oi and 02 parameters suggested by the data, remained a theoretical 
puzzle. (Even before the advent of QCD factorization various authors had presented 
arguments against the universality of nonfactorizable effects in class-2 decays; see, e.g., 
||13| , |T4| , |T5|].) New experimental data announced by the CLEO and Belle Collaborations 



T6|, |T7| change the picture significantly, in a way that is entirely in line with QCD 



expectations. As we will illustrate below, these data provide compelling evidence for 
process-dependent 02 values, and for large relative strong-interaction phases between 
related class-1 and class-2 amplitudes. This shows that the "naive" factorization model 
is too simple to account for the data. 

Because QCD factorization cannot be justified for the class-2 decays, it is in some 
sense misleading to parameterize the — >■ D'^ir^ decay amplitude as done in (|l|). Al- 
though the naive factorization contribution 02 = C2 + (Ci (with ( a nonperturbative 
parameter of order 1/Nc) is certainly present, it is not a leading contribution to the decay 
amplitude in any consistent limit of QCD. For instance, there exists a weak annihilation 
contribution to 02 which scales like a power of rub/A in the heavy-quark limit and thus 
formally dominates over the naive factorization piece. (Weak annihilation is mentioned 
here only as an example of a leading contribution to the class-2 amplitude. Model cal- 
culations suggest that the annihilation contribution in B ^ Dn decays is nevertheless 
small 0. Another example of a leading contribution is charge-exchanging rescattering 
from the dominant class-1 channel WR, 0.) Also, for class-2 decays naive factorization 
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Table 1: Heavy-quark scaling laws in the schemes where the charm quark is 
treated as a heavy quark (left) or as a light quark (right) 





rric ^ rrib ^ rriQ rric ~ A ^ 


h/2 — h/2 


l + 0(A/mQ) 0(1) 
0(A/mQ) 0[asimh)] 



does not emerge in the large- iVc limit. It is then more appropriate to employ an alterna- 
tive parameterization of the decay amplitudes in terms of isospin amplitudes A1/2 and 
corresponding to transitions into Dtt final states with / = | and |, respectively. It 
is given by 




AiB- ^ D\-) = Vs A3/2 . (3) 

An identical decomposition holds for other decays such as -B ^ D*n and B D'^*^p. It 
follows from QCD factorization that the ratio of isospin amplitudes is 

"^'/^ -l + 0(A/mQ), (4) 



V2Ay2 

which also implies that the relative strong-interaction phase 61/2 — ^3/2 = 0{A/mQ). 
Here mq represents either one of rric and m^. Note that the corrections to the "1" in (^) 
are also formally suppressed by a power of l/N^, and it has been argued that perhaps 
this color suppression may be more relevant to factorization than the heavy-quark limit 
| 20| . However, an identical 1/Nc argument would apply to other nonleptonic decays such 



as D ^ K^*^7r and K ^ rni, for which color suppression is clearly not operative. Apart 
from trivial substitutions of quark flavors, the only difference between, say, B D^*^n 
and D K^*^7r decays is the larger energy release in the decay of a heavy b quark, which 
leads to color transparency and thus is the basis of QCD factorization. 

The deviation of the ratio A1/2/ {V2 A3/2) from 1 is a measure of the departure from 
the heavy-quark limit. When contemplating about the expected magnitude of this effect, 
it is important to realize that the power suppression of the corrections in (^) relies on 
the heaviness of the charm quark, not only the b quark. In order to illustrate this fact 
it is instructive to consider two different power-counting schemes for the heavy-quark 
expansion. The most natural scheme, which underlies our discussion so far, is to consider 
both beauty and charm as heavy quarks with their mass ratio rric/ f^b fixed in the heavy- 
quark limit. Then the deviation of Ai/2/("\/2 ^3/2) from 1 is power suppressed in A/mg, 
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Table 2: Experimental data for the B D^*^7i and D K^*^7!' branching 
ratios (in units of 10~^), isospin amplitudes, and related quantities 





B^Dn B ^ D*7i 




D ^ K-K D ^ K*7f 


B- dWOtt- 


3.0 ±0.4 2.76 ±0.21 
0.27 ±0.05 0.17 ±0.05 
5.3 ±0.5 4.6 ±0.4 


DO K^*)-n+ 


38.3 ±0.9 50 ±4 
21.1±2.1 31±4 
28.9 ±2.6 19.0 ±1.9 


IA1/2/V2A3/2I 

^1/2 — ^3/2 


0.70 ±0.11 0.72 ±0.08 
(27 ±7)° (21 ±8)° 


IA1/2/V2A3/2I 

51/2 — (^3/2 


2.71 ±0.14 3.97 ±0.25 
(90 ±6)° (104 ±13)° 


X 02/01 
arg(a2/ai) 
xaf/af 

af/af 


0.42 ±0.05 0.35 ±0.05 
(56 ±20)° (51 ±20)° 
0.25 ±0.12 0.23 ±0.08 
^0.28 ^0.25 


X 02/01 
arg(a2/ai) 
xaf/af 
af/af 


1.05 ±0.05 1.11 ±0.08 
(149 ±2)° (160 ±2)° 
-0.53 ±0.02 -0.66 ±0.02 
~ -0.44 ^ -0.35 



where niQ ~ mj, ~ m^. It is not calculable without model dependence. Alternatively, 
one may consider the charm quark as a light quark with its mass kept fixed in the heavy- 
quark limit IQ. In such a scheme the class-2 amplitude becomes calculable and of leading 
power in the limit m^, ^ A . The leading deviation of the isospin amplitude ratio from 
1 is then computable in terms of Wilson coefficient functions Ci{mb) and short-distance 
corrections proportional to a<j('^fe)- The scaling of the relevant quantities in these two 
versions of the heavy-quark limit is summarized in Table The conventional scheme is, 
perhaps, somewhat closer to reality. However, considering that the charm quark is not 
very heavy in the real world, we may expect a sizeable deviation of the amplitude ratio 
from 1. 

We now investigate to what extent the prediction (|) is supported by the data on i? ^ 
0^*^71 decays. For comparison, it will be instructive to analyze the related charm decays 
TT in parallel. The upper portion in Table summarizes the experimental data 
on the various branching ratios. The color-suppressed B decays have just been observed 
for the first time experimentally |16|, |1^]. The preliminary results for the B^ D^tt^ 
branching ratio are (2.6 ± 0.3 ± 0.6) x 10"^ (CLEO) and (2.9 Ig;^ ± 0.6) x 10"^ (Belle), 
while those for the 5° D*°tt° branching ratio are (2.0 ±0.5 ±0.7) x 10"^ (CLEO) and 
(1.51o;51q;4) x 10~^ (Belle). We have averaged these results to obtain the entries shown 
in the table. All other numbers are taken from P2|. By combining the measurements 



of the three branching ratios for each mode, taking into account the lifetime ratios 
t{B-)/t{B^) = 1.068 ± 0.016 and t{D+)/t{D^) = 2.547 ± 0.036 pi, one can 



extract the magnitude and phase of the ratio of isospin amplitudes. The results are 
shown in the middle portion of the table. A clear qualitative difference between beauty 
and charm decays emerges. Whereas the relative phases are close to maximal in D decays 
and the amplitude ratios are far from the asymptotic value in (H), the corrections to the 
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Figure 1: Allowed regions in the complex y4i/2/(A/2 ^3/2) plane obtained 
at 68%, 95% and 99% confidence level. The sign of the imaginary part is 
undetermined by the data. The black dots show the central values. 



heavy-quark limit appear to be much smaller in B decays. 

A more careful analysis of the experimental data for B ^ Dn and D Kir decays is 
shown in Figure |l], which gives the allowed regions for the ratio v4i/2/(v^ A3/2) obtained 
at different confidence levels. To derive these regions we find, for each value of the isospin 
amplitude ratio, the minimum of the function for the three measured branching ratios. 
We then plot contours of minimum iii the complex plane corresponding to a given 
confidence level (for two degrees of freedom) . Very similar constraints can be derived for 
B D*TT and D —* K*tt decays. Note that even at the level of one standard deviation 
(68% confidence level) the relative strong phase of the isospin amplitudes for B decays 
shown in the left plot may be zero, in contrast with the naive error propagation in 
Table 0. In the right plot we compare the results in D and B decays and also indicate 
the value corresponding to the strict heavy-quark limit (HQL). We conclude that the 
data on the ratio of isospin amplitudes in B D^*\ decays is compatible with the 
heavy-quark scaling laws discussed earlier. The comparison of charm to beauty decays 
shows a clear progression towards the heavy-quark limit as the mass of the decaying 
quark increases. The remaining deviation of the amplitude ratio from 1 is compatible 
with a correction whose suppression is governed by the large charm-quark mass. 

We have mentioned earlier that the parameterization of the class-2 amplitude in 
(H) is somewhat misleading, because the naive factorization contribution is in no way 
the leading term in a controlled expansion of the decay amplitude. Nevertheless, it 
is instructive to extract the parameter a2 defined via the second relation in (|l]) and 
compare it with the values of the 02 parameters obtained from other decays, such as 
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Figure 2: Allowed regions in the complex xaijdx plane for B 0^*^71 decays, 
obtained at 68%, 95% and 99% confidence level. The sign of the imaginary 
part is undetermined by the data. The black dots show the central values. 



B K^*\j/ip. The ratios of the B D^*^it branching ratios can be expressed in terms 
of the ratios xa2/ai, where 



x(Dtt) = ^—^ ■^^ " ^09 



The numerical values have been obtained using f£, ^ 200 MeV, f^* ~ 230 MeV, and 
F(f^'^(m|))/F(f ^^(m^) ^ F^^'^ [m]^,) / A^^^* {m\) 0.5. (The corresponding quan- 
tities for D decays may be estimated using the BSW model [|I|, with the result that 
x{Ktt) 1.2 and x{K*Ti) ^ 1.9.) Figure |] shows the corresponding allowed regions 
for xa^jax obtained at different confidence levels. The central values are also shown in 
Table 0. The data prefer values \a2\ ~ 0.4-0.5, which are larger by almost a factor 2 than 
those obtained from B — > K^*^J/il) decays (see above). [| Moreover, the best fits prefer 
phases of about 50° in the two cases (with large errors), suggesting that strong final- 
state interactions cannot be neglected and lead to a nontrivial relative phase between 
class-1 and class-2 amplitudes. Both observation are in conflict with the assumptions 
underlying the naive factorization model. It now appears that the 02 coefficients of the 
class-2 amplitudes are nonuniversal, with magnitudes and phases that may be rather 

^Preliminary Belle data p7[ | on related color-suppressed decays support this conclusion. Using the 
form- factor model of HI, we find that the decays B° D^*^°-q and B° D'-*'^"uj have \a2\ values 
between 0.4 and 0.5, with experimental errors of about 0.1 and a theoretical uncertainty of about 30%. 
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different for different types of decays. More precise data will be required to fully explore 
the pattern of QCD effects in the class-2 and class-3 amplitudes. 

We should mention that also in D — ^ K^*"*!: decays the extracted values of 02/01 have 
large phases and are larger in magnitude than those extracted from other D decays (see 
Table This "failure" of naive factorization is usually attributed to the strong final- 
state interactions caused by nearby resonances (as signaled by the fact that |5i/2 — 53/2I ~ 
90° in these decays). It is then argued that one can only expect to correctly predict the 
magnitudes of the isospin amplitudes but not their relative phase [Q] . The ratio of these 
magnitudes is determined by the ratio of the real, effective a^^ and a'^ parameters of 
the naive factorization model via the relation 



IA/2I _ 2-xaf/af 



721^3/2! 2{l + xaf/af) 

The effective parameters so determined are shown in the lower portion of Table In the 
case of D decays, the physical picture underlying this approach is that of predominantly 
elastic final-state interactions, which mix the various K^*"*!! final states and thereby 
changes the phases but not the magnitudes of the isospin amplitudes. While this as- 
sumption may be questioned even in the case of charm decays |2^, is it clearly not jus- 



tifiable for decays of B mesons, in which rescattering is predominantly inelastic ^ . 
Therefore, we believe it is a coincidence that the "effective" a"^ /a^ ratios are close to 
the expectations of the naive factorization model. 

Finally, we like to comment on the observation that the ratios a; in (^, which govern 
the relative strength of the class-2 and class-1 amplitudes in naive factorization, exhibit 
large violations of the scaling x ~ (A/mQ)'^ expected in the heavy-quark limit. Does this 
imply a failure of QCD factorization in hadronic B decays? We believe the answer to 
this question is negative. Consider first the conventional case where the charm quark is 
treated as a heavy quark. Then the ratios x arise only in naive factorization. The fact 
that they are not numerically suppressed reflects the well-known failure of the conven- 
tional heavy-quark expansion for heavy- light form factors and decay constants, i.e., the 
empirical fact that the ratios 

^-1.5 h(A/mQ)V^], h(A/mQp] (7) 

do not scale as expected from the heavy-quark limit of QCD. (The reason for this failure 
may be related to the "smallness" of /^r, which in turn reflects the smallness of the light 
quark masses via the relation = —2{mu + m^) {qq).) However, as we have argued 

earlier the factorized class-2 contributions appearing in the numerator of the ratios x 
are not a leading contribution to the class-2 amplitudes in the heavy-quark limit. Other 
contributions exist that are parametrically larger. It is therefore not clear to what 
extent the large scaling violations in (0) are relevant to the class-2 amplitudes. In the 
opposite limit where the charm quark is considered a light quark the naive factorization 
contribution is the leading contribution to the class-2 amplitudes. In this limit also the 
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ratios x are of leading order in the heavy-quark expansion, which is consistent with the 
numerical values x{D^*'>t^) ^ 0.9. So there is no evidence for a failure of the heavy-quark 
expansion either. 

In summary, we have argued that new experimental results on the color-suppressed 
nonleptonic decays — > £)(*)°7r° provide evidence for a failure of the naive factorization 
model and for sizeable relative strong-interaction phases between class-1 and class-2 
B — > D^*^TT decay amplitudes. This resolves a long-standing puzzle created by the 
apparent universality and small rescattering phases of the class-2 parameters 02 in B 
decays. The new data suggest that the 02 parameters in different types of decays such 
as S — > D^*^ TT and B — > K'^*^J/ijj differ by almost a factor 2 in magnitude, indicating a 
strong nonuniversality of nonfactorizable effects. This is in agreement with theoretical 
expectations based on the heavy-quark expansion. We find that the size of corrections 
to the heavy-quark limit seen in the data is compatible with the expectation that the 
suppression of the corrections is governed by the large charm-quark mass. We urge our 
experimental colleagues to produce more precise data on a large variety of hadronic B 
decays. This will help to explore in detail the pattern of QCD effects in the class-2 and 
class-3 amplitudes, as well as to further establish the validity of QCD factorization (and 
hence the applicability of the heavy-quark limit) for class-1 decays. 

Acknowledgement: This work was supported in part by the National Science Foundation. 

Note added: While this Letter was in writing the paper hep-ph/0107257 by Z.-z. Xing 
appeared, in which a similar analysis is carried out but different conclusions are obtained. 
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